Temperature dependence of the plasmonic Casimir interaction 
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We investigate the role of surface plasmons in the electromagnetic Casimir effect at finite temperature, in- 
cluding situations out of global thermal equilibrium. The free energy is calculated analytically and expanded for 
different regimes of distances and temperatures. Similar to the zero-temperature case, the interaction changes 
from attraction to repulsion with distance. Thermal effects are shown to be negligible for small plate separations 
and at room temperature, but become dominant and repulsive at large values of these parameters. In configura- 
tions out of global thermal equilibrium, we show that the selective excitation of surface plasmons can create a 
repulsive Casimir force between metal plates. 

PACS numbers: 73.20.Mf - Surface plasmons, 05.30.-d - Quantum statistical mechanics, 68.35.Md Surface thermodynamics, 
surface energies in surfaces and interfaces 
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I. INTRODUCTION 

The interaction between two parallel plates due to the zero 
point fluctuations of the electromagnetic field is commonly 
known as the Casimir effect. For metallic plates, it is well 
known that at short distance, the interaction can be attributed 
to surface plasmon modes 111 01 that hybridize across the vac- 
uum between the interfaces JslQl- Surface plasmons have 
been attracting much interest in the last years in connec- 
tion with a broad range of topics such as near-field spec- 
troscopy, sub-wavelength resolution Q |g] or extraordinary 
optical transmission through subwavelength metallic hole ar- 
rays I7H9l|. The electromagnetic field associated with these 
modes is evanescent. It therefore came as a surprise that they 
also give a large contribution at large distances, and even a re- 
pulsive one Ill0l - [l2ll . These papers have been restricted to the 
Casimir effect at zero temperature. The present paper gen- 
eralizes these results by including a nonzero temperature and 
situations out of thermal equilibrium. One might expect that 
the thermal excitation of surface plasmons is irrelevant, since 
their typical energies are comparable to the plasma frequency 
of the metal, much larger than experimentally relevant ther- 
mal energies. As retardation is taken into account, however, 
the surface plasmon dispersion relation approaches the light 
cone and drops to lower frequencies. These are comparable, 
for two parallel plates, to the lowest cavity resonance ~ c/L. 
We find indeed a significant thermal component to the Casimir 
interaction between surface plasmon modes when the distance 
L exceeds the thermal wavelength ^ hc/ksT. By selectively 
exciting a class of plasmonic modes, we even get an overall 
repulsive Casimir force. 

The paper is organized as follows. We first recall the dis- 
persion relations for coupled surface plasmon modes on two 
metallic plates llll and obtain a general expression for the cor- 
responding Casimir free energy. This is expanded asymptot- 
ically in different regimes of distance L and temperature T 
in Sec. [in] Sec. |IV] discusses the plasmonic Casimir entropy. 
We then compare the results to the full Casimir interaction 



between metal plates (Sec.|V]i, including all electromagnetic 
modes, consider situations in which the plasmonic modes are 
not at the same temperature as the rest of the system (SeclVIli, 
and conclude with a short summary. 

We adopt throughout this paper the following lossless di- 
electric function 
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where cUp is the plasma frequency. The Casimir energy and 
force (both per unit area) are normalized to the values found 
for perfectly reflecting mirrors 11311 
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where H = ISO/tt'^. An intrinsic physical length scale of the 
system is the plasma wavelength Ap = 2TTc/ujp. It is con- 
venient to use the latter as a length scale, switching to a di- 
mensionless plate distance A = L/\p. We introduce also a 
reduced temperature r = T/Tp = {2TTkBT)/hwp, where Tp 
is the plasma temperature. This choice makes the numeri- 
cal results independent of the specific material, and gives uni- 
versal scaling laws. Note that 1/r is proportional to the ra- 
tio between the thermal wavelength (a few microns at room 
temperature) and the plasma wavelength, and that the product 
At = kBTL/{hc) is independent of the plasma wavelength. 
A parameter set used frequently in related work is the one for 
gold: hup = 8.96 cV, Ap == 136 nm and Tp = 1.66 x 10*K. 
Room temperature then corresponds torwl.SxlO^^. 



II. PLASMONIC CASIMIR FREE ENERGY 

The Casimir free energy of two metallic plates is obtained 
by summing the free energy (per mode) over the electromag- 
netic modes vibrating inside the cavity yllTsI]. This expres- 
sion is suitably regularized, namely by subtracting the limit of 
large distances between the plates. 




Figure 1 : Dispersion relation of the surface plasmon modes at plate 
separation L = c/uj^ (dotted curves) and L = 2c/ajp (dashed 
curves); k is the wavevector parallel to the surface. The solid black 
curve u}o{k) corresponds to the surface plasmon on an isolated plate 
(L — >■ oo). The branch (^+ (k) crosses the light cone. The asymptotic 
value at large k is uj^p — ujf/y/2. 



Since the modes of the electromagnetic field are formally 
equivalent to harmonic oscillators, the free energy of a sin- 
gle mode of frequency ui in thermodynamical equilibrium at 
temperature T is 



free energies of the coupled modes and subtracting twice the 
corresponding values at infinite distance the integral over the 
dispersion relations gives the plasmonic Casimir free energy 
in the form (coefficients c± — \, cq = —2) 
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which is convergent at large k llOtl . The thermal part of 
the free energy [second term in Eq. (|3]l] naturally cuts off 
modes above ksT /h. Different choices of the lower limits 
ka, related to the subtraction procedure, are possible and have 

been discussed in Refs. iflU [l2l [l4l - [l6ll . They are connected 
with the way the evanescent and propagating contributions of 
the mode ui+ {k) are split. Here we apply the convention of 
Refs. Illll[l2ll and set ka = for all modes. We thus include 
both propagating and evanescent branches of the 'plasmonic 
mode' w-|_(fc). 

The calculation of the integral (01 is challenging because 
the surface plasmon dispersion relations uta {k) are solutions 
of a transcendental equation, except in the non-retarded limit 
k ^ ojp/c where 0Ja{k) ~ Wsp(l + ae '^^ with ojsp = 
ujp/\/2, a ~ 0,±. Progress can be made with the paramet- 
ric form described in Refs. Ill0l - [l2ll . Adopting the notation of 
Ref. lUin . we get the dispersion relations LOa{k), a e {0, ±}, 
from 



^l{z) = ji^gliz), 
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Kiz) = j^ (5) 
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In this paper, we sum Eq. (O over the dispersion relations 
for the surface plasmone modes. Isolated surfaces (at infi- 
nite distance) carry a single surface plasmon mode of fre- 
quency ajo(fc), illustrated in Fig. [T] If the plates are brought to- 
gether, the electromagnetic fields of the modes overlap, break- 
ing the degeneracy and splitting the dispersion relation in two 
branches, whose frequencies we label Lu±{k). The modes 
uj^{k) and ujo{k) are both entirely evanescent and lie be- 
low the light cone. The mode aj+(fc), however, crosses the 
light cone and connects smoothly with the lowest propagating 
mode (with p-polarization) within the cavity [1]. Adding the 



with the dimensionless functions 
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(The exponents are ±1 for a = ±.) The parameter z varies 
from —Za . . . oo: one has zq = z^ = 0, while the number z+, 
plotted in Fig.|2]as a function of A, is the solution 1 1111 of the 
transcendental equation ^/z^ = 27rAcos(yzIj7/2). The prop- 
agating branch of the mode uj+ (k) corresponds to the interval 
z G [— z+, 0]; evanescent modes (below the light cone) have 
z > 0. Changing the integration variable in Eq. dUi from k to 
z, the plasmonic Casimir free energy is given by 



J 



T 



he 

87rL3 



Yl ""^ 



a=±fi 



oo 

/ (gaiz) 



2Ar In 



1 — e >-^ 



9a (oo) 



dz + 2 



{gl + ^9a 



Arln 



1 - e ^^ 



dga 



9a(-7-.) 



(7) 



For all plasmonic modes, the second integral has a finite up- per limit (7a(oo) = v^ttA = Wgpi/c that coincides with 



the non-retarded surface plasmon frequency; their contribu- 
tions cancel in the sum over a. The lower boundaries are 
9o{-zo) = 9-{-z_) = 0, and.g+(-z+) = ^/zj. The first 
integral in Eq. (|7]i can only be evaluated approximately (see 
Sec.lnili. but a closed form can be given for the second one. 

In the following, we scale the plasmonic free energy to the 
zero-temperature Casimir value, Eq. (|2]i. 



T{L,T) = Ec{L)^{\t) 
and split the correction factor in two terms 

^(A,T)=r?(A)+^(A,r), 
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where the first is the plasmonic Casimir energy at zero tem- 
perature ifTHfTl : 



VW = -2^^^" 9a{z)dz + 



3/2 



(10a) 



The second term in Eq. (|9]l gives the temperature-dependent 
part for which Eq.© gives 



^(A,T)=-HAT^Ca 
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where the following combination of polylogarithmic func- 
tions appears 



£(x)-C(3)-Li3(e-")-xLi2(e-^) 



(11) 



with Li„(.T) = J2T=i x''/k". We note that C{x) - ^x^{l - 
21oga;) for small x, and C{x) — )■ C(3) exponentially fast for 
large x. Eq. (llObb does not depend only on the product At be- 
cause the material-dependent parameter A enters via the lower 
limit z+ and the functions ga{z) [Eq. (|6]l]. 




Figure 2: 2+ (A) vs. the plate separation A and its asymptotes at 
small and large distances. The limiting cases are 2+ ~ (27rA)'^ and 
« -K^ for A — >■ and -^ 00, respectively. 



0.13 




/ 


/^ 


-^ 




0.10 
0.05 


/ 




y 


\^ 


2\ ^ 

i.sX \ - 


0.00 


J'^ 






\ 


\ \ \ 








\ 


\'\ \ \ 


0.05 








r 


^ \\\\ ^ 


o.in 










\ \ \\ ; 



Figure 3: Plasmonic contribution to the Casimir free energy vs. dis- 
tance at different temperatures, normalized to the perfect mirror case 
at T = (energy correction factor v^(A, r) in Eq.([8}). Solid curves: 
numerical evaluation of Eqs. jlOal llObt for different temperatures; 
dashed curves: short-distance limit Eqs. l |15l|16t . Distance and tem- 
perature are scaled to the plasma wavelength 2-nc/u]^ and tempera- 
ture Tp = hwp/kB, respectively. Negative values correspond to a 
repulsive interaction energy. 
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Figure 4: Thermal correction to the plasmonic free energy at short 
distance for different temperatures (energy correction factor ^^(A, r) 
in Eq.([9l(). The behavior qualitatively differs from the complete free 
energy (Fig.[3}, as a quadratic (rather than linear) distance depen- 
dence emerges at low temperatures, cf. Eq. J16t . 



The Casimir free energy is shown in Figs. [3] IH 111 and [8] for 
different distance ranges. Qualitatively, a nonzero tempera- 
ture does not modify the behavior of the plasmonic contribu- 
tion - we still get a sign change at a distance of order Ap/27r, 
with the interaction becoming repulsive at large distances. We 
plot in Fig. |5] the inversion distance where the Casimir pres- 
sure, —dF/dL, changes sign: a weak increase is found as the 
temperature is raised. Much larger changes will be found in 
Sec|Vl] where configurations out of thermal equilibrium are 
discussed. 

In the following, we analyze the thermal correction ^(A, r) 
in different regimes of distance and temperature. The zero- 
temperature Casimir energy 77(A) depends only on one physi- 
cal scale provided by the plasma wavelength and leads to two 
regimes A ^ 1 and A 3> 1. For ??(A,r) we discuss three 
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Figure 5: Repulsive and attractive regimes of the plasmonic Casimir 
pressure in the (A, r)-plane. A nonzero temperature slightly in- 
creases the distance for which the plasmonic contribution to the 
Casimir force becomes repulsive. 



regimes. In all situations of practical interest, it is safe to 
assume r ^ 1 for the scaled temperature and we can distin- 
guish between: short distances A ^ 1, intermediate distances, 
1 ^ A ^ 1/r, and large distances (beyond the thermal wave- 
length) 1/r < A. 



III. ASYMPTOTIC EXPRESSIONS 



A. Short distance 

At short distance the zero-temperature energy correction 
was already analyzed in Ref. Ill2ll . It turns out that it is dom- 
inated by large values of z (third integral of Eq. (fT4]i). At the 
leading order, we get 



^(A) 
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> 1.790A. 



(15) 



Higher order terms take the form A'' (a + & log A) with numer- 
ical coefficients a and h given in the same reference. 

Considering the thermal correction for A <C 1, the ther- 
mal scale becomes important, too. For realistic temperatures, 
we also have At <C 1, and the main contribution arises from 
the second and the third integral in Eq. ( fT4l i. Indeed, it can 
be shown that the first integral and the polylogarithmic term 
(involving £) are beyond the order 0(}?\ In addition, the 
main contribution to the second integral arises from the mode 
u!-{k). This is not surprising since the thermal correction se- 
lects frequencies uja{k) < T and the mode uj-{k) is the one 
that vanishes most quickly as fc -^ 0. The corresponding ex- 
ponent in k determines the power law in r, as we discuss at 
the end of this section. The opposite case Ar ^ 1 is phys- 
ically irrelevant at short distances, because one would need 
T ^ 1. Mathematically, one finds a divergence from the term 
proportional to C in Eq. (llObb that is exactly balanced by the 
first integral in Eq. (fT4l i. It follows that the asymptotic form 
given in Eq. (fT6] l remains valid. 

All told, up to the second order in A we find 



The calculation of asymptotic expressions requires some 
care already at zero temperature as was shown in Ref. uM- 
When performing approximations on the integrals in 
Eqs. dlOab and (llObb one must bear in mind that the func- 
tions ga{z) cover a wide range of values, from very small 
to large, depending on z and A. Their characteristic scale 
in the variable z is given by the distance parameter A. For 
— Za < 2; <C (27rA)^ we may use (a = 0, ±) 
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while for z 3> (27rA)^ the (non-retarded) approximations 
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hold. It is therefore convenient to split the integration range in 
Eqs. ( llOab and ( llObb as follows 
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The first integral concerns only the mode cj_|-(fc) because 
z_ , zo = 0. We can use Eq. (fTSl i in the first two integrals 
and Eq. ( fTST l in the third. Depending on distance and temper- 
ature and on the desired accuracy, these have to be compared 
with the integrated terms in Eq.dTol) (proportional to z^ or 
/:(V^/Ar)). 
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where the function C{x) defined in Eq. (fTTT i appears with a dif- 
ferent argument, and where the temperature-dependent func- 
tion /3(r) is 
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This is plotted in Fig.|6] together with its asymptotes in the 
limits of high and low temperatures. 
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We observe the emergence of the characteristic ratio A/r^ 
that determines which of the two terms in Eq. ( fT6b dominates. 
This illustrates that the limits A ^- and r — )• do not 
commute for the temperature-dependent Casimir energy. If 
A <C r^ (extremely short distances or high temperatures, main 
plot of Fig.|4|i, the function /3{t) governs the thermal correc- 
tion -d which scales as -d ^ Xt^ if A ^ t^ ^ 1. We recover 
here the same linear distance dependence as at zero temper- 
ature, Eq. ( fTSl ). The opposite regime r'^ <C A ^ 1 emerges 



B. Intermediate distance 



at low temperatures, where the term involving L in Eq. (fTSI l 
dominates: we get a behavior -d ~ }?t^ (inset of Fig.|4]i. This 
crossover from a quadratic to a linear scaling with distance 
can be seen in Fig.lD 



This discussion also 
retarded approximation, 
persion relations uj\ (fc) : 



illustrates the failure of the non- 
This leads to surface plasmon dis- 

and 



(1 ± e~^^) and wq 



"sp 



extrapolates a free energy •& ^ t' down to low temperatures, 
while the correct power is t^ . This is of course crucial for the 
low-temperature expansion of a thermodynamic quantity like 
the entropy (see Section UVb . 
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Figure 6: The function /?(r) vs. r and its asymptotes ( 118b at low and 
high temperatures. 



Parameters for typical experiments are r « 10^^, A « 
1 . . . 10^, way beyond short distances. They lie inside an in- 
termediate regime 1 <C A ^ 1/t, where the plate distance 
is between the plasma and the thermal wavelength. Here, the 
thermal correction to the free energy is still small compared 
to zero temperature, as for short distances. 

In the scaled Casimir energy r\{\^ r), the main contribution 
to the integral ( I lOab arises for z ^ \ and we have to con- 
sider the first two integrals in Eq. (fT4l) . where, approximately, 
z^ « Tv^. The energy correction factor at zero temperature 



then becomes ll 
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A>1 



-74.57\/A + 60. 
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the offset arising from the second term in Eq. dlOab . 

For the thermal correction, we find the leading order from 
the polylogarithmic L in Eq. jlObb . The z-integral gives 
a contribution which is dominated by the interval z = 
— z+...(27rA)^. Combining the two. 
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This gives a small correction that scales as t^ . Both Eqs. (fT9l l 
and (|20] | are plotted in Fig.|7] illustrating the weak impact of 
temperature. It is interesting to note that in this range of dis- 
tances, the thermal plasmonic contribution is opposite in sign 
to the free energy of the full electromagnetic Casimir effect 
(see Sec. lIVI ). and increases the plasmonic repulsion. 



C. Large distance 
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Figure 7; Thermal plasmonic Casimir free energy (reduction factor 
<^(A, r) in Eq.©) at intermediate distances and different tempera- 
tures. Exact numerical calculation (solid curves) and the approx- 
imation ( 119b at zero temperature (short-dashed white) and sum of 
Eqs.(|19b and ( 120b in the intermediate regime \t <^\ (long-dashed) 
respectively . 



Let us finally consider the regime A ^ 1/r ^ 1, cor- 
responding to a plate separation larger than both the plasma 
and the thermal wavelength. The zero-temperature contribu- 
tion can still be approximated by Eq. ( fT9] l, but now the thermal 
contribution dominates the free Casimir energy. The asymp- 
totic behavior of the integrals in Eq. dlObb is obtained by ex- 
panding the logarithms for small (7(j(z)/(At), since the func- 
tions 5a (^) are bounded: 
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Using this expansion under the integral in Eq. ( llQbb . we 
note that the second term balances exactly with the zero- 
temperature contribution from Eq. (llOal ). As for the first term, 
we perform the z-integration by splitting the integration range 
as in Eq. (fl4b . It is easy to see that in the second interval, the 
sum over the mode branches gives zero. The leading contri- 
bution now comes from negative z, while the third interval in 
Eq. ( fT4] i gives an exponentially small contribution ^ er^^^. 
The polylogarithmic term L can be expanded for small argu- 




only the branch of the plasmonic mode w+(fc) that crosses 
into the propagating sector contributes to the (repulsive) plas- 
monic Casimir interaction at large distances in a significant 
way. 



IV. PLASMONIC CASIMIR ENTROPY 

The plasmonic Casimir entropy can be derived from the 
plasmonic Casimir free energy, Eq. (|7]i, by differentiation with 
respect to T, 



Figure 8: Plasmonic Casimir free energy (correction factor ip{X, r)) 
at large distances and different temperatures. Numerical calculation 
(solid curves) and large-distance approximations J19! ) at zero temper- 
ature (white dotted line) and {22} at non-zero temperature (dashed 
curves). 
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A convenient scale is given by the Casimir entropy at high 
temperatures between two perfect reflectors (this includes two 
transverse photon polarizations) 
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where the dependence on In t cancels to leading order The 
validity range of this asymptotic formula is illustrated in Fig. [8] 
(dashed lines) where the full free energy is plotted at large 
distances. 

Summarizing this section, we have gen eralized a result 
known from the zero-temperature case lUll [l2r to T > 0: 
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C(3) kB 
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The scaled entropy is connected to the dimensionless thermal 
correction t9(A, r) by a derivative (mind that Eq < 0) 
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We recall the result for perfect reflectors where the entropy 
depends only on the product Ar iflTIl 



'7c(A,r) 
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1, At>1. 



(26) 



The Casimir entropy due to surface plasmons can be repre- 
sented as the integral 



a{X,T) 
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where na{z) = [exp{ga{z) / \t) — 1]^^ is the Bose-Einstein 
mean photon number. Figs.|9] [To]show the temperature de- 
pendence of a{\, t) for several distances A, small and large. 
The strong qualitative differences between these cases could 
be anticipated from the free energy of Fig.lH at short dis- 
tances, temperature makes increase i?(A,r) towards positive 
values, leading to a positive a from Eq. dZST l, while the trend 
is reversed at larger distances. 

Note from Figs|9] and [TO] that the plasmonic Casimir en- 
tropy fulfills Nernst's heat theorem (ct ^- as T — > 0) at all 



distances. Though the same is known for the entropy of the 
complete plasma model lUSll including photonic modes, the 
result is not trivial because surface plasmons are only a sub- 
system of the two-plate system. 

As before one can distinguish three characteristic distance 
regimes for the plasmonic Casimir entropy. The expression 
for A <C 1 can be easily obtained from the approximation (fT6b 
to i?(A, t). At low temperatures (regime r^ ^ A ^ 1), we 
must include a subleading term in the small-A expansion to 
get the right prefactor of the temperature power law. This is 
done by adding to Eq. ( fTSI l the poly logarithmic term with C of 
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The validity range of these formulas can be seen from Fig.fTOl 
(dashed lines). Note that for intermediate and large distances, 
the entropy approaches zero from below as r — > 0. 



Figure 9: Temperature dependence of the Casimir entropy from plas- 
monic modes [correction factor cr(A, r) relative to perfectly conduct- 
ing mirrors, Eq. ( I23H for short distances and high-temperature limit 
J28b (dashed). Inset (double logarithmic scale): low-temperature be- 
havior at short distance and its asymptote (dashed) cr ~ r^ from 
Eq.l[27ll. 
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Figure 10: Plasmonic Casimir entropy in the scaled form u{\, r) 
[Eq. (I23H. vs. temperature for intermediate and large distances (solid 
lines). Dashed: low- and high-temperature asymptotes J29b . J30b . 



Eq. ( llQbb that becomes - 
to 

a(r,A)^ac(A,r) 



-2H(At)3C(3). Differentiation leads 
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n^ C(3) 



T 



(27) 



The entropy approaches zero quadratically as r ^ 0, as for 
perfect reflectors [Eq. (|26]l1. but the prefactor is larger by a 
factor l/(7rA). The good agreement with the exact result can 
be seen in the inset of Fig.|9] 

At high temperature, ^(A, t) becomes a constant that coin- 
cides at short distances with the perfect reflector limit for one 
polarization 



a(A, 






(28) 



This can be seen from Eq. ( fTSI l taking into account the func- 
tion /3(t). 

Intermediate (1 < A < 1/t) and large (1, 1/r < A) dis- 
tances can be treated with Eqs. (|20] |. ( |22] | that give, respec- 



V. PLASMONIC VS. PHOTONIC MODES 

We now compare the plasmonic Casimir interaction to the 
full system where all electromagnetic modes of the cavity are 
included. The knowledge of the surface plasmon contribu- 
tion provides some physical interpretation for the complete 
system. For example, it is well known ||3l 01 that the full 
Casimir interaction (zero temperature) at small plate separa- 
tions is well described by taking only the electrostatic inter- 
action between surface plasmons. We now show that this re- 
mains valid at nonzero temperature. The fundamental reason 
is that at short distances (A ^ 1), the lowest cavity modes 
(above the plasmonic uj^ + (fc)) have a characteristic fre- 
quency 2ttc/L that already falls in the transparency band of 
the mirrors (ui > Wp). 

The full Casimir free energy J^if = V'Lif (A, t)Ec can be 
obtained from the Lifshitz formula lUSl [1911 . In our scaled 
units. 



<^Lif(A,T) = -2HAT^^'rp(27rnAT), (31) 

p n—0 

TpiX) = / dKKlog[l-r2(iX,K)e-2«]. (32) 
Jx 

The index p G {TE, TM} denotes the polarization. The num- 
bers Xn = 27rnAT are scaled Matsubara frequencies. The 
Fresnel reflection coefficients in terms of the variables k and 
X are 



rTEiT^X,K) 



rTui'iX^K) 



e{iX)K — K„ 
e{iX)K, 



(33) 

with K,n = a/k^ + (27rA)2 and the dielectric function of the 
plasma model [cf. Eq. ^] e{iX) = 1 + (27rA/X)2. 

Fig.[TT] shows the scaled free energies for both the full 
Casimir interaction (dashed) and the plasmonic contribution 
alone (solid). A good agreement is visible at short distances 
even at nonzero temperature. In this regime, we can, there- 
fore, obtain detailed information on the thermodynamics of 
the Casimir effect by just considering the plasmonic contribu- 
tion which can be worked out in analytic form quite easily. For 
example, we can immediately conclude that the (full) Casimir 
entropy at A ^ 1 is given by Eq. (|27| i at low temperatures 
(r < 1) and by Eq. (I28]l for t > 1. 
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Figure 1 1 : Casimir free interaction energy including all modes (Lif- 
shitz theory with the plasma model, dashed lines), compared to its 
plasmonic counterpart (solid). All energies are expressed relative via 
the correction factor f{\, r). 
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Figure 12: Casimir entropy vs. temperature for all modes (plasma 
model, dashed), compared to the contribution of plasmonic modes 
only (solid). Entropies expressed via the correction factor a(\, r). 



With respect to the Casimir entropy. Figs. [TTI and [TSlillus- 
trate that the Lifshitz expression deviates significantly from 
the plasmonic contribution when A > 0.1. It is obvious that 
propagating (photonic) modes then become relevant. In Ap- 
pendixlAl we calculate their contribution at low temperatures 
for the plasma model: 



T<1 



Stt^t 2 + ttX' 
45C(3) 3 



(34) 



where the first term is twice the value obtained for perfect 
mirrors. This is precisely compensated by the plasmonic con- 
tribution. Indeed, for intermediate distances 1 <C A <C 1/r, 
both Eqs. (|29] l and ( l34b are valid, and their sum reproduces the 
entropy of the full Casimir effect calculated in Refs. lUSl |20| - 
I22II . Evaluating Eq. (19) of Ref. 1I21I1 in the regime of interme- 
diate distances, we have 



cTLiflA, r) -^ — >ac{X,T) 



1 



SttV ttX + 2 



45C(3) 



VI. BEYOND THERMAL EQUILIBRIUM 

Until now, we have assumed both metallic slabs to be at the 
same temperature T. The previous results enable us to deal in 
a simple way with a more general situation, too, where each 
of the (otherwise identical) slabs is described by a local tem- 
perature Ti and 72. The general theory in this case was inves- 
tigated in Refs. ll23ll24ll : the non-equilibrium Casimir interac- 
tion for a symmetric cavity is obtained by simply averaging 
over the equilibrium free energies of the two mirrors 



.F'-q(L,T2,ri) = i[^<=q(L,T2) 



.F^'^(i,Ti)] 



(36) 



(35) 



(In Ref. |24f| , this result was derived for the pressure, but the 
same reasoning can be applied for the free energy.) Combin- 
ing this formalism with the results from the present paper, it is 
straightforward to calculate the plasmonic contribution to the 
non-equilibrium Casimir free energy. Using the split ^ of the 
plasmonic free energy into a zero-temperature and a thermal 
part, Eq. ( l36l l gives for two slabs at different temperatures 

^"^'^(A,r2,Ti) - ?/(A) + i [^(A,r2) +^(A,Ti)] . (37) 

From the results given above, we conclude that qualitatively, 
the behavior of iy9"'"^(A, T2, ti) is similar to the equilibrium 
configuration, including a change in sign of the force with the 
distance. This is also confirmed by the asymptotic expressions 
for long/short distance and low/high temperature that can be 
easily extracted from the above results. The total Casimir 
force between identical plates, however, is always attractive, 
as is known from Refs. ll23ll24ll for all temperatures. 

Let us now consider a slightly different non-equilibrium 
scenario where temperature is still raised locally (in one 
plate), but only for a subclass of modes. If it were possi- 
ble to increase the mean excitation of the plasmonic modes 
on one plate, above the equilibrium level of the propagating 
(photonic) modes, the total Casimir free energy would read 

¥'u'(A,r2,ri) = ^uf{X,n) + i [z?(A,T2) - d{X,n)] , 

(38) 
where the first term is the total Casimir free energy at equi- 
librium. The set of curves b) in Fig. [13] illustrates that this 
scenario can create a regime where the total Casimir force 
becomes repulsive, and this over a fairly large range of dis- 
tances. We plot the Casimir pressure (non-equilibrium force 
per unit area) when the photonic modes are either at the scaled 
temperature ti (zero or room temperature), and the plasmonic 
modes on plate 2 at T2 > ti . It appears that this setting breaks 
the delicate balance between photonic and plasmonic modes 
we found in Sec.|V] A similar interpretation has been put for- 
ward in Ref. 12511 for the change in distance dependence of the 
atom-surface interaction out of equilibrium. The two values 
for Ti give close results because in the intermediate distance 
range, the effect of the temperature is still moderate for the 
equilibrium case. As could be expected, the inversion dis- 
tance increases and the maximal repulsion becomes weaker 
as n increases towards T2. 
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The selective excitation of surface plasmon modes is a well- 
studied problem (Ref. 12611 and references therein). Most of 
the setups have to cope with the fact that the corresponding 
electromagnetic field is evanescent and, therefore, cannot be 
excited directly by laser photons incident from free space. 
Corrugated surfaces are of some help here 131, and indeed they 
can convert thermally excited plasmons into far-field radiation 
luTIl . A recently developed four-wave mixing scherne permits 
to excite surface plasmons even on flat surfaces ll26ll . The 
non-equilibrium situation involving only propagating modes 
[Eq. dTTT il may be simpler to realize experimentally since these 
modes couple to free-space light fields and can in principle be 
excited by shining a laser I28|| from the side onto the gap be- 
tween the mirrors. (See Ref. 112911 for a related discussion.) 



Figure 13: Total Casimir force (per unit area) in thermal equilibrium 
(thick line) and in different non-equilibrium scenarios. The force 
is normalized to 10~®|-Fc(Ap)|, approximately 3.65 /iPa for gold 
(Ap = 136nm). In these units, T = 300 K (665 K) corresponds to 
r Ri 0.018 (0.04), respectively. 

a) Total equilibrium at temperature r: attractive pressure at all dis- 
tances, b) Surface plasmon modes of one plate out of equilibrium at 
temperature r2, all other modes at temperature t\. c) All modes at 
r, except for the propagating branch of the plasmonic mode (aj+ (fc)) 
which is at temperature Tp,-. The sign change to repulsion (positive 
pressure) would occur for gold at distances between ~ 2.7 /xm and 
3.7/im. 



We have also included in Fig.[T3]a scenario where only the 
propagating part of the plasmonic mode w+(fc) is populated 
at a temperature different from the rest of the system. It con- 
tributes a free energy 



^P'(A,r)=-H(Ar) 



111 



2H {XtY 



C 



Xt 



(39) 



L 



.g+(0) 
At 



where 5+ (- 



z+ as noted before, and 



5+(0) = 2^A 



1 



I + ttA 



(40) 



gives the (dimensionless) wavevector for which the dispersion 
relation cj+(fc) crosses the light cone. This leads to a non- 
equilibrium free energy 

C'(A,V,r) =^Lif(A,T)+z9P'-(A,rp,)-z9P!(A,T) , (41) 

where Tpi is the temperature of the propagating plasmons and 
r is the temperature of all other modes. The corresponding 
pressure [curve c) in Fig. [131 increases with respect to the pre- 
vious non-equilibrium scenario by approximately a factor of 
2, and repulsion sets in at a somewhat shorter distance. This is 
because (i) the - otherwise attractive - mode w_ (fc) is less ex- 
cited and (ii) the propagating branch of a;+(fc) dominates the 
interaction at these distances and excites the electron plasma 
on both plates rather than a single one. 



VII. DISCUSSION AND CONCLUSION 

We have calculated the contribution to the thermal Casimir 
effect due to surface plasmons, which are hybrid field-matter 
eigenmodes of metallic surfaces. The expression we found for 
the free energy of interaction is valid at any distance and tem- 
perature, and we have derived its asymptotics at small, inter- 
mediate, and large distances. Thermal effects become signif- 
icant when the distance is larger than the thermal wavelength 
At, similar to perfectly conducting plates, and below A^ for 
non-equilibrium configurations. The other length scale of the 
system (plasma wavelength) determines the detailed behavior 
of the free energy. 

We have found that at short distances and temperatures the 
thermal coiTection is small, and that the plasmonic Casimir 
interaction changes sign with distance, leading to a repulsive 
regime, as has been known from zero temperature lUU [VM . 
This goes hand in hand with a change of sign of the plasmonic 
Casimir entropy. In the short-distance regime, we found that 
the complete Casimir interaction between metallic plates (de- 
scribed by the plasma dielectric function) is completely dom- 
inated by the surface plasmon contribution. The asymptotic 
scaling laws explain why T = is a good approximation in 
most experimentally relevant situations (intermediate distance 
regime, low temperature). In this regime the known result for 
the complete plasma model is recovered in a simple way by 
adding propagative photonic modes. 

Things are different at high temperatures and large dis- 
tances. Here it was shown that the plasmonic Casimir interac- 
tion is determined by a branch of the surface plasmon disper- 
sion relation corresponding to propagating modes, resulting in 
a large repulsive contribution that is enhanced by the temper- 
ature. This effect is probably one of the best illustrations of 
Casimir repulsion that arises from the radiation pressure of a 
standing wave mode. The pressure is repulsive, because the 
travelling photons are bouncing off the cavity walls, while the 
reference mode, a single-interface plasmon, has an evanescent 
field with zero radiation pressure. 

The balance between plasmonic and photonic modes was 
emphasized by considering two configurations out of global 
thermal equilibrium where plasmonic modes are selectively 
excited to a higher temperature. These configurations show 
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a crossover to a total Casimir force that becomes repulsive 
at plate distances L w 20 Ap . . . 25 Ap (a few microns for 
gold). This can be understood qualitatively in terms of radia- 
tion pressure due to the propagating branch of the plasmonic 
mode. We emphasize that this happens at distances shorter 
than the thermal wavelength where the Casimir pressure is 
stronger 

In conclusion, it seems in principle possible to tune the sign 
of the Casimir force by the selective excitation of the surface 
plasmons. Still, future research must address experimentally 
relevant questions for such a scheme, e.g. how to avoid excit- 
ing photonic modes just above the plasmonic one and how to 
populate plasmonic modes over a wide angular range. 

Acknowledgments. We would like to thank H. T. Dinani 
and S. Slama for discussions and help with some calculations. 
We benefited from exchanging ideas within the Research Net- 
work "Casimir" of the European Science Foundation (ESF). 
FI acknowledges partial financial support by the Humboldt 
foundation and LANL. HH and CH acknowledge funding by 
the German-Israeli Foundation for Development and Research 
(GIF). 



The first (second) integral in Eq. (IA2b corresponds to the 
evanescent wave (propagating wave) sector, respectively. For 
p = TE, the argument of the logarithm is always positive in 
the first integral, hence its imaginary part vanishes. This does 
not happen for p = TM where the first integral gives the con- 
tribution of surface plasmons (evanescent branch) which has 
been evaluated in this paper. As mentioned in Sec.|Vl we are 
interested here in the propagating contribution only. 

The function x/ sinh x significantly differs from zero only 
for X < 1. In the limit Xt ^ 1, we can therefore expand the 
integrands in A'Ip{n) for small y and 51 since y < H, -^ I. 
This yields 



4 OTTA 



ImM^l^in) 



SttA 



-n^ + ^^ {3 + nX) 



(A3) 
(A4) 



Performing the a--integral in Eq. ( lAll i. 



Appendix A: Full Casimir entropy at low temperatures 

The dimensionless correction factor for the Casimir entropy 
of the plasma model can be written as the following integral 
over (scaled) real frequencies 

4 / X dx X — ^ 

f7Lif(A,r)= — - / _-^ImVMp(2a:;AT) (Al) 

I'U'^J Jo smh X ^ 

where p = TE, TM indicates again the polarization and 

/•OO 

Mp{n)^ dKKln[l-rl{n,K)e-^''] 

Jo 

+ / dyy\n[l-rl{n,-iy)e^'y] (A2) 

^0 



(Tph(A, t) -^ — > a-c(A, r) 

T<1 



Stt^t 2 + 7rA 
45C(3) 3 



(A5) 



Note, however, that this result contains the propagating 
branch of the plasmonic mode w+(fc) whose free energy is 
given by Eq. (|39] l. Reviewing the analysis from Sees. IIII Al 
nil Bl it is easy to see that in the limit considered here, 
the polylogarithmic term dominates in Eq. ( [39] l and becomes 
£(y'zIj7/Ar) « C(3)- Subtracting this contribution from ( IA5I ). 
we find the entropy of the propagating photonic modes given 
in Eq. dUl. 
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